function I = simpson( a, b, n, f )
%	Functia care calculeaza valoarea integralei functiei f pe
%intervalul [a; b] , integrala calculata prin metoda Simpson.
%
%	Date de intrare:
%		- a -> extremitatea stanga a intervalului( pe care se calculeaza integrala);
%		- b -> extremitatea dreapta a intervalului( pe care se calculeaza integrala);
%		- n -> gradul de precizie considerat;
%		- f -> functia a carei integrala se cere de calculat.
%
%	Date de iesire:
%		- I -> valoarea integralei.

h =(b-a) / (2*n);
s1 = 0;
s2 = 0;
for i = 1 : n
	s1 = s1 + f(a+(2*i-1)*h);
endfor
for i = 1 : n-1
	s2 = s2 + f(a+2*i*h);
endfor
I = h*(f(a) + f(b) + 4*s1 + 2*s2)/3;

endfunction
